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An unified thermodynamical framework based in the use of a generalized Massieu-Planck thermo- 
dynamic potential is proposed and a new formulation of Boltzmann-Gibbs Statistical Mechanics is 
established. Under this philosophy a generalization of (classical) Boltzmann-Gibbs thermostatistics 
is suggested and connected to recent nonextensive statistics formulations. This is accomplished by 
defining a convenient squeezing function which restricts among the collections of Boltzmann-Gibbs 
configurations of the complete equilibrium closure. The formalism embodies Beck-Cohen super- 
statistics and a direct connection with the nonlinear kinetic theory due to Kaniadakis is provided, 
being the treatment presented fully consistent with it. As an example Tsallis nonextensive statistics 
is completely rebuilt into our formulation adding new insights (zeroth law of thermodynamics, non 
ad hoc definition of the mean value of a physical quantity,...). We relate all the formal development 
to physical and measurable quantities and suggest a way to establish the relevant statistics of any 
system based on determinations of temperature. 
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I. INTRODUCTION 

The physical nature of theprobabilistic approach in Equilibrium Statistical Mechanics has been discussed exten- 
sively in the literature (see P, S IS S IE and references therein) . The entropy has long been considered as the 
relevant function connecting probability theory with the thermodynamical description of the physical world. Jaynes 
proposed in 1957 a formulation of Statistical Mechanics based on the so-called Maximum Statistical Entropy Prin- 
ciple within the framework of Information Theory, by identifying the thermodynamic entropy S with the maximum 
value of the statistical entropy functional C{{pi}) {pi denoting the probability of an event i) through the Boltzmann 
constant k. This is accomplished by means of the relation S/k — C. Interesting criticisms of Jaynes approach have 
been clearly pointed out in the references (see Q , and references therein) . 

In the last two decades, there have appeared several probabilistic views with their grounds on Jaynes approach, 
but applied to new definitions of entropy. These methods have their origin in the pioneering work by Tsallis 
who unfolded the theory of nonextensive statistics of great subsequent interest. This theory has been succesfuUy 
applied to a wide variety of systems and some points of the formalism are currently still under debate. Besides 
Ruelle's previous works (for a recent and clear discussion see, for example, and references therein) and Hill's 
nanothermodynamics [Tol lllj| , nonextensive thermostatistics contain new interesting views on systems that go beyond 
traditional thermostatistics incorporating fractality, complexity and strong correlations. The existence of these systems 
are nor surprising nor exotic, and has originated a debate on the validity of the ergodic hypothesis which has been, 
on the other hand, reinforced and extended after rigorous works by other authors pj- There is no reason, however, 
to believe a priori that all systems in nature are ergodic. Boltzmann, in fact, encountered examples of clearly 
non-ergodic systems but considered them as exceptional, unimportant, cases (the well known example of the linear 
chain of coupled oscillators being an example). The study of nonequilibrium situations and complex systems has 
thus motivated the apparition of new formalisms that constitute important efforts to face the inherent complexity of 
long range interactions, long term memory, fractality, etc. There exists however some open questions and difficulties 
concerning nonextensive thermostatistics as for example the definition of thermal equilibrium. 

In this paper we propose a new approach to both equilibrium and nonequilibrium (stationary states) statistical 
mechanics. We also address the problem of providing a sound thermodynamic basis to the statistical methods 
presented which has its roots in a generalized version of the zeroth law. Applications of the theory to concrete 
physical systems and formal developments will be given in forthcoming papers but we advance that the methods 
proposed can be applied to a wide variety of systems provided that concrete developments as Tsallis statistics, Beck- 
Cohen superstatistics and Hill's nanothermodynamics are embodied in the present formulation. We begin discussing 
equilibrium statistical mechanics and, in the last section, we take advantage of the philosophy developed at equilibrium 
to extend the formalism to nonequilibrium situations and complex systems. Our postulates are no longer explicitly 
referred to the entropy but to a more general thermodynamic function which attains a minimum at equilibrium: the 
dimensionless characteristic thermodynamic function A", which coincides with —S/k for an isolated system and with 
the Helmholtz free energy in kT units F/kT for a system closed with diathermal walls in contact with a heat bath at 
temperature T. X is always dependent on environment variables through the generalized partition function which is 
regarded as the collection of possible (weighted) configurations compatible with the constraints. After introducing two 
postulates, a theorem concerning the form of X is proven. We have, for any system 

X = -\ng (1) 

where g is the generalized partition function of the system. We shall no longer use the concept of partition function, 
however, and we consider g as the collection of configurations within a class specified by a choice of independent 
variables (see below) . If the system is isolated then the collection of configurations g coincides with the microcanonical 
(countable) configurations. In this case Eq.QJ takes the form of the celebrated expression for the microcanonical 
entropy due to Boltzmann, S = fclnfi, where f2 are the microcanonical configurations (i.e. the total number of 
attainable microstates) of the system. 

The condition of minimum for the appropriate thermodynamic potential for a system at equilibrium is discussed 
extensively in the literature (see, for example. Landau and Lifshitz' book T^l, §20 and §21). It implies that the 
minimum work required to drive a system from an equilibrium state to a neighboring one is always positive. We make 
use of this point. This work is defined as the increment of the thermodynamic characteristic function (= kTAX) when 
intensive variables are kept constant. The application of the condition for a minimum in X leads us to the central 
expressions of ensemble theory. If all attainable configurations within the classes are allowed, Boltzmann-Gibbs (BG) 
equilibrium Statistical Mechanics is derived. We consider that X depends on the variables chosen to characterize the 
system (which we shall call environment variables) through the number of attainable configurations by the system. 

In a recent paper by Vives and Planes jl3j |. the importance of Massieu-Planck thermodynamic potentials in con- 
necting Thermodynamics with Statistical Mechanics is emphasized. Our work shares some connections with some 
points of the work by these authors, but the approach is quite different because the authors essentially follow the 



3 



maximum entropy principle as formulated from information theory. An advantage of our formulation, besides the 
compactness achieved in the expressions (formally identical for all physical systems whose microcanonical classes (see 
below) are known) is that it provides an interface between thermodynamics and statistical mechanics that can be 
easily extended to non-conventional situations as self-organized systems or systems capable of taking decisions. In the 
last section we propose a scheme for the generalization of BG equilibrium statistics. As a straightforward example, 
Tsallis non-extensive statistics is reconsidered and rebuilt within this scheme by introducing only one additional 
assumption. 

The formulation presented here is general and applies also to small systems (with the only exception of Section 
III El which deals with fluctuations in macroscopic systems) when homogeneity of X is not demanded. We summarize 
below some features and issues of the formulation that could be of general interest: 

1. Compactness and generality in the expressions derived. Only the knowledge of the dimensionless characteristic 
thermodynamic function X in terms of generalized extensive variables Xj and intensive ones yi are required to provide 
a complete thermodynamic description of any system. If the configurations g are also considered in terms of these 
variables an statistical framework can be developed. 

2. Hill's nanothermodynamics 0,^3 is embodied in the approach followed. This formalism (developed previously) 
is, besides Gross microcanonical thermostatistics 14] (not discussed here), a quite convenient framework to explain 
the nonextensive (equilibrium) thermodynamics of systems with size comparable to the range of the forces involved. 
It is important to make here a distinction between an extensive or nonextensive thermodynamics and an extensive 
or nonextensive statistics. A nonextensive statistics does not necessarily imply a nonextensive thermodynamics. Im- 
portantly, Hill's formalism makes stable the previous classical development of Thermodynamics to the domain where 
fluctuations in physical quantities are important. 

4. The statistical formulation can be easily extended and applied to non conventional situations (complex and 
self-organized systems, nonequilibrium and metaequilibrium, etc.) The method is general and it is not restricted, in 
principle, to some concrete choices made below to illustrate how it works. 

5. In particular, Tsallis statistical mechanics is rebuilt within our framework. New insights concerning this theory 
are provided. Concretely, a new definition for the mean value of a mechanical quantity is suggested. This definition 
is different to those encountered in the foundational papers of nonextensive statistics (see below) and is not given ad 
hoc but arises naturally from the formalism. Tsallis form for the entropy in terms of probabilities is derived. 

5. The generalized formulation presented can also be related easily to recent developments on Statistical Mechanics. 
We discuss, concretely, Beck-Cohen superstatistics. 

6. The formulation is fully consistent with the kinetic theory previously proposed by Kaniadakis which generalizes 
the Boltzmann kinetic theory and the Fokker-Planck picture. Besides this previous work, irreversibility can be 
understood in a broader context. The close interrelation between both formalisms makes the formulation presented 
well grounded kinetically and can thus serve as a basis for the study of a wide variety of statistical systems. 

The paper is organized as follows: in Section III Al the dimensionless characteristic thermodynamic function X is 
introduced from thermodynamical considerations showing how the knowledge of this function provides the complete 
thermodynamical description of the system. In Section fll BI two postulates on the dimensionless characteristic function 
are presented, relating it to the configurations that a physical system can attain. An important theorem on the form 
of the function X is also proven. In Section III CI ensemble theory is developed over the two preceding postulates. 
In Section III Dl probabilitv is introduced and the Gibbs-Shannon form for the entropy is derived. In Section III El a 
general theory of fluctuations is given. In Section III Fl the previous formalism is applied to the usual ensembles. 
Finally, in Section IlIII it is shown in general how BG equilibrium statistics can be extended to complex systems and 
nonequilibrium. Squeezed statistical mechanics is formulated and a connection with the nonlinear kinetic theory by 
Kaniadakis, which provides a generalization of the celebrated Boltzmann kinetic equation, is established. In Section 
IIVI Tsallis nonextensive statistics is discussed and completely rebuilt as a concrete formal application by introducing 
an appropiate squeezing function motivated from physical grounds. New insights in Tsallis statistics are presented. 
The connection with Beck-Cohen superstatistics is also established. 

We remark that, although Sections 111 Al to Hi El are largely conventional (as are referred to classical BG equilibrium 
thermostatistics) the constructive procedure that we use to develop the non conventional statistics (Section}^) is 
presented. 
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II. FORMULATION OF BG THERMOSTATISTICS 



A. The characteristic thermodynamic function as a thermodynamic potential 

For any given system, the following differential form for the entropy J — S/k holds 

dJ~Yl yi'^^i = (2) 
I 

Here Xi is an extensive variable and yi is its conjugate intensive one divided by kT. For a simple pure substance under 
hydrostatic conditions {Xi} = {£',¥, iV} and {yi] = {l/kT,p/kT, —fi/kT} are the two sets of variables and Eq.Q 
(where E, V, N, p and fi are the energy, volume, number of particles, pressure and chemical potential, respectively) 
can be written as 

In Eq.(|3J the differential of the entropy is expressed in terms of the differentials of E, V and N being this the adequate 
form (for practical purposes) to describe the thermodynamics of a completely closed system (as it is characterized by 
the specification of these three extensive variables). However, for a partially open system, other choices of variables 
are more adequate. The variables chosen to characterize the system are called environment (or natural) variables. 

Let us consider the sets of variables {Xi} and {yi}- Each of them can be separated into two subsets satisfying 
{Xi} = {Xj}U{Xi} and {yi} = {yj}^{yi}- We define {Xj} as the subset of extensive environment variables and {yi} 
as the subset of intensive environment variables. The subsets {Xi} and {yj} contain the conjugate variables (which 
are not environment variables) of those within the former subsets. Thus, an environment extensive variable Xj is the 
conjugate of an intensive non-environment one yj, while an intensive environment variable yi is the conjugate of an 
extensive non-environment one Xi. Hence, the complete set of environment variables is {Xj} U {yi}. For an isolated 
system {yi} = {0} and, therefore, {Xi} = {0}. 

For a specification of M environment variables, if the subset {Xj} contains m variables, {yi} contains M — m 
variables and Eq. (j^J can be rewritten as 

m AI 



dJ-Y.y^'^^i^ J2 y^dX.^O (4) 

j — 1 i—m+1 

If we add and substract a term d (^YldLm+i Ui-^ij > obtain 

(J\/ \ m M 

J- y^^^ ) - yj'^^i + = ^ (5) 

t— m+l / J — 1 i—m+1 

We define the dimensionless characteristic thermodynamic function X as 

X= y,X,-J (6) 

i—m+l 

We adopt henceforth Einstein summation convention for repeated indexes (and then X = yiXi — J). From Eq.jS)) we 
thus have 

dX = -yjdXj + X,dy^ (7) 

This differential equation can be considered as a generalized form of Eq. Q which was written for a completely closed 
system. As we have seen above, for an isolated system {Xj} — {Xi}, {yj} — {j/;}, {Xi} — {0}, {yi} = {0}, and, 
from Eq. ©, X — —J. With these considerations, Eq.jSJ is recovered. It is important emphasizing here that 
Eq.0 is valid for any system independently of its constraints. This equation also shows an important feature of the 
dimensionless characteristic function: given the set of environment variables, the knowledge of X in terms of these 
allows the complete thermodynamical description of the system. Non-environment variables yk {yk G {Vj}) ^nd X^ 
[Xs e {Xi}) are then specified through the following relations 



Xs ^ { — ] (9) 



dysj 
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The dimensionless characteristic function X is a generahzed Massieu-Planck function. Massieu-Planck functions 
are entropic thermodynamic potentials defined as Legendre transformations of the entropy |13| . For the concrete 
cases of an isothermal closed system and an isothermal-isobaric impermeable system, the dimensionless characteristic 
function defined by Eq.® coincides with the thermodynamic potentials first introduced by Massieu and Planck 
[Id] , respectively through a minus sign. 

In macroscopic systems X is an homogeneous function of first order in extensive environment variables and of 
zero order in intensive environment variables. Hence, due to the linearity of the differential form of X, Eq. Q, the 
following relation holds 

X = -y,X, (10) 

for any macroscopic system. Ea. HlU|) is the well known Euler's equation for a macroscopic system. As a consequence of 
it, by taking the differential in each side of Ea. (|l()(l . and comparing the result with Eq.Q we obtain the Gibbs-Duhem 
equation for macroscopic systems 

Xjdy^+X,dy,^0 (11) 

One can relax, in general, the requirement of X to be an homogeneous function of first order in extensive environment 
variables. This is the case of small systems where this homogeneity breaks. An outstanding feature of these systems 
is the presence of finite-size effects due to forces of range comparable to the size of the system and to its coupling to 
the environment. Hill's nanothermodynamics |lOl lllj is an important approach to explain their thermal properties. 
Here we follow a different approach which proves to be equivalent to Hill considerations. We suggest that for small 
systems the deviations from the behavior of macroscopic systems in the thermodynamic limit can be accounted for 
by introducing a term O in Ea. HlUf) by the following way 

X = -yjXj - O (12) 

The Gibbs-Duhem equation takes now the form 

dO = -X,dy, - X,dy, (13) 

These considerations are consistent with Hill's definition of the subdivision potential £ [lol[ll| . In fact O = ~E/kT. 
Note that within this way of introducing the inhomogeneities, formal symmetries appear in the equations of thermo- 
dynamics among variables chosen to characterize the system. On one hand Eqs.Q and (jl3|l 

dj = yidXi (14) 
dO = -Xidyi (15) 



and on the other Eqs.® and 



X = y,X, -J = -y,X, - O (16) 



show that there exists two different entropic functions O and that can lead, through the convenient Legendre 
transforms, to the same thermodynamic description of a system, specified by the dimensionless characteristic ther- 
modynamic function X . The dimensionless entropy corresponds to a thermodynamic description in which the 
independent variables are all extensive. In the case of a description based on O, the independent variables are all 
intensive. For a partially open system, the appropriate function is X . The conventional entropy of the system is ob- 
tained in terms of X by closing completely the system through the adequate Legendre transforms while O is obtained 
opening completely the system by making the Legendre transforms in the reverse order. 

The results of this section have shown that the knowledge of X in terms of the environment variables is enough 
to completely describe the thermodynamical behavior of a system. In the following section X becomes the central 
function in connecting the microscopic "counting over configurations" to the macroscopic world. It is important 
to remark at this point that thermodynamical non environment variables constitute observed values of microscopic 
averages. Establishing the form that these microscopic averages take is one of the objectives of the following sections. 
If we denote by obs thermodynamically observed variables, Eq. H17|l can be rewritten as 

X = y^X.^obs - J (17) 
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B. Postulates of Equilibrium Statistical Mechanics 

A physical system is found, in a given instant of time, to be in a particular configuration. Constraints imposed 
to the system limit the number of configurations that it can adopt. When some of these constraints are removed 
the system evolves through an irreversible process occupying new possible configurations until equilibrium (where 
new but less restrictive constraints are specified) is reached. The result is that the total number of configurations 
has been increased attaining a maximum on going to equilibrium, where the function X attains a minimum. These 
facts always hold for an irreversible process. A larger number of attainable configurations is related to a lower value 
of X . Furthermore, each configuration contains a partial thermodynamic description in itself: it is compatible with 
the specification of some variables. We can then introduce the concept of collections or classes of configurations 
9{x...}u{y...} to denote the sets of configurations compatible with the specification of a set of variables. On collecting 
these classes of configurations we obtain a complete thermodynamic description of the system if we relate them to 
a thermodynamic potential. The dimensionless characteristic thermodynamic function X is the most general form 
for an entropic thermodynamic potential, embodying in itself all the choices of concrete entropic thermodynamic 
potentials that can be done, compatible with specifying a set of environment variables. 

For an irreversible process at constant intensive environment variables we have, from Ea. (|17|) 

AX = y,AX., -AJ<0 (18) 

for the system on going to equilibrium. Eq. (|18|l implies that it is necessary to do a positive work on a system to drive 
it from equilibrium to a nonequilibrium neighboring state [T^ . This consideration applies, in principle, to systems 
with many particles, as it is of thermodynamic character. It is not obvious to hold for systems with few particles due 
to the microscopic reversibility of the dynamical behaviour [T^. However, it constitutes a good starting point for all 
what follows since it assures the very desirable property of stability (convexity in the free energy, concavity in the 
entropy) that can be extended later to other situations. This demanded convexity in the free energy is not a matter 
for worrying about when studying, for example, phase transitions of small systems with the apparition of convex 
intruders in the entropy, since we supplement the information concerning small systems (see below) throu gh t he Hill's 
subdivision thermodynamic potential (the theory for this important concept can be found elsewhere [lOlllll|'). 

It is useful to consider in more detail the concept of classes of configurations. At equilibrium these classes can 
be hierarchically interrelated and defined, the classes being unions of subclasses. In general, the most important 
class is what we call below the characteristic class of the system: the collection of elements compatible with the 
specification of environment variables {Xj} and {yi}, g{Xj}u{yi}- Another important class is that of the total number 
of microstates of the objective system, the microcanonical class g{x.j} iii which only the extensive variables of the 
system are specified. For an isolated system, the microcanonical class is also the characteristic one. In equilibrium, 
all the classes and interrelations between them can be established, and the main classes described above can be 
decomposed in subclasses through the following way 

9{x,} = U 9{x,}u{x,} = 9{x,}u{x,} (19) 
9{x,}u{v,} = U 9{x,}u{v,}u{x,} = 9{x,}u{y,}u{x,} (20) 

{X,} {Xi} 

where the sums are over all the possible specifications of each of the variables belonging to the set {Xi}. The 
microcanonical subclasses 9{Xj}u{Xi} are particularly interesting since they constitute the so-called degeneracy of a 
configuration, and can be evaluated directly from the Hamiltonian after specifying variables {Xj} and {A^i}- Knowing 
how the characteristic class is decomposed in microcanonical subclasses is, hence, mandatory and, stated in other way, 
it constitutes the celebrated Gibbs theorem. We present below a new derivation of this theorem that has a certain 
(vague) resemblance with the foundations of statistical mechanics by the method of counting in the microcanonical 
ensemble due to Balian and Balazs [3- We choose, however, to focus the discussion in the concept of classes of 
configurations, in which a system can be, instead of probability. This is appropiate when studying the important issue 
of irreversibility since microcanonical classes evolve with time, a feature which is absent in probability distributions 

An irreversible process can be pictured as the result of removing one or more constraints. Removing these implies 
enlarging the collection of configurations as the classes of configurations are formed by elements compatible with the 
constraints and for less restrictive constraints there exists more elements compatible with them. It is important to 
note here that this statement holds for systems which are uncapable of restricting (or enlarging) by themselves their 
collections of configurations. The role of the constraints is crucial (and the only essential thing) for these systems in 
constructing the classes. In the last section of this paper we study complex, organized systems that can select among 
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restricted (or enlarged) classes of configurations. Let us study first the case of systems that have not this capability, 
i.e. Boltzmann-Gibbs (and, as we shall see, ergodic) systems. We postulate: 

Postulate 1 For any physical system there exists a dimensionless characteristic thermodynamic potential X which 
is a monotonously decreasing function of the size of the characteristic collection of configurations of the system 
(abbreviated by g) i.e. X(g) < X{g+ \ Ag\). g coincides with the numerical value of the appropiate partition function. 
Therefore, when the system is isolated g coincides with the total number of microstates (the microcanonical class). 

Postulate 2 The dimensionless characteristic thermodynamic potential is an extensive quantity, i.e. for two inde- 
pendent systems A and B, Xa+b = Xa + Xb- 

These two postulates allow us to build the whole formalism of ensemble theory and fluctuations. This is shown in 
the following sections. We prove first two theorems on the form of the dimensionless characteristic thermodynamic 
function. 

Theorem 1 For independent identical systems, each of them with g elements X satisfies X{g-^) = NX{g). 

The proof of this theorem is immediate if one considers that, if the systems are independent, the number of elements 
of the whole ensemble are given by g^. Let be the dimensionless characteristic function of the ensemble and Xi 
that of one system. Then, by Postulate El we have Xj^ = -^(5^) = A/'A'i = NX{g). 

Theorem 2 The dimensionless characteristic thermodynamic function is given by X{g) — —\\\\D.g where A is a 
dimensionless constant. 

In the proof of this theorem we follow an analogous treatment to the developed by Landsberg for the entropy 
(pp. 126-128 of his book). Let r, n be arbitrary positive integers, with m determined by the following chain of 
inequalities 

< g" < r™+^ (21) 

Forming \nx for each term of this equation, one finds 

rn In r < n\ng < (m + 1) In r (22) 

whence 

m In g m -|- 1 , , 

- < < (23) 

n inr n 

An exactly similar relation is satisfied by A", by virtue of Theorem ^ 

m X(n) m + 1 , , 



It follows that 



Xig) In 5 



X{r) Inr 

The l.h.s. of Ea. (|25|l is independent of m, while n is arbitrary. Hence, upon taking n arbitrarily large, it is obtained 
that 

m _ m (26) 



< - (25) 
n 



\ng In r 
is a constant. A, independent of r. Then 

X{g)^\\ng (27) 

From Postulate^we see that A < and, hence, we obtain the desired result. 

|A| is a dimensionless constant which adjust the scale of temperatures. For example, for a mean number < > of 
molecules of ideal gas occupying a macroscopic volume V of an isothermal open region one finds X = —pV/kT = A < 
iV >. The equation of state of the ideal gases can be reproduced empirically and the fitting of the experimental curves 
leads to an expression pV =< N > kT for the ideal gases. Then, A can be conveniently identified with a numerical 
value of —1 and, hence, we can write 

X{g)^-\ng (28) 

We shall use this result in all what follows. 
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C. Ensemble Theory 

The three fundamental expressions obtained for X, Eqs. (|17|l . Q and (|28() contain, expressed in compact form, the 
main expressions of formal ensemble theory of equilibrium statistical mechanics and thermodynamics. A workable 
version of the former in the present framework is necessary and this is the aim of this Section. Ea. H28|l can be 
rewritten, in terms of the specified environment variables, as 

X = -\ng[x,}u{y,} (29) 

The main advantage of focusing the discussion in the function X is that all equations obtained are valid for any system. 
Let us consider a completely isolated system. All environment variables for this system are extensive and these specify 
a total number of configurations 5{Xj} = 9{Xi} = ^ compatible with the constraints imposed. For this system the 
dimensionless characteristic thermodynamic function X coincides with — J', see Ea. (jl7|l . The (dimensionless) entropy 
then can be defined as minus the microcanonical dimensionless characteristic function. This definition of the entropy 
is general and valid also for partially opened systems if one always considers the whole number of configurations 
compatible with regarding the system as completely closed. As for a partially opened system there are different ways 
of closing it all configurations compatible with the different ways of closing the system should be counted and added. 
In the following we consider a discrete spectrum of configurations for simplicity. Analytically we can write for the 
entropy J in general 

^ = In X! 9{x,}u{x,} = lng{x,} (30) 

where the property stablished in Eq. (|19|l has been used. Note that this definition is referred to the total number of 
configurations compatible with extensive environment variables. The very meaning of closing a system is specifying a 
concrete value for the extensive variables. The set of extensive environment variables change for different systems, but 
the definition of the entropy is always formally identical. A constraint is imposed, in general, when a thermodynamic 
variable is forced to have a concrete value. As indicated above, microcanonical classes are directly related to the 
Hamiltonian and, hence, it is convenient to know the combination of classes of configurations g{Xj}u{Xi} on which 
9{x }u{yi} is decomposed. 

We have, from Eqs.^ and (gnj) 

X ^- \ng^x,}u{y,} = - In XI 9{x,}u{y,}u{x,} (31) 

If we specify a concrete value for the extensive non environment variables we have, by using Eqs. (|17|l and H3U|) 

■^{x,} = Vi^i - J{Xi} = ViXi - h^g{x,}u{Xi} 

= -ln[5{x,}u{x.}e~^'^'] (32) 

where Xi_obs — Xi has been used in Ea. H17|) since this variable is now specified. From Ea. (|31|l . we can write 

^{Xi} = -lng{Xj}u{yi}u{Xi} (33) 
and, by comparing Eqs.|j22Jl and H33|) we obtain 

9{x,}u{y,}u{x,} = g{A',}u{x,}e"^'^' (34) 

and therefore, from Eq.lPU 

X=-\nY,9{x,}u{x^}e'y'''' (35) 

Hence, Gibbs theorem is proven and ensemble theory is thus formulated through the previous approach. The arguments 
presented above have lead us to general expressions obeyed for all systems in Boltzmann-Gibbs equilibrium. We 
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summarize below the most important equations obtained: 

?i = -ln.9{^^.}u{y,} 

dX = -yj^obsdXj + Xi^obsdjji (37) 

J = —X + HiXi^obs 

= ln.9{Xj} (38) 
J-{xj - -X[x,^+y,X, (39) 

The subindex i* makes reference to aU variables belonging to the set {yi} excluding only that involved in the derivative. 
Note that in Ea. (|36|l the sum is carried over the whole spectrum of values that extensive non- environment variables 
can take, as these are allowed to vary freely. Thus, Ea. (|36|l serves as a definition of X. The argument in the logarithm 
is the generalized partition function. Ea. H37|l allows the thermodynamical description of the system and Ea. (|38|l is a 
definition of entropy as a function of X. Finally Ea. (|39|l is a useful definition of the entropy when concrete classes of 
configurations compatible with the specification of non-environment extensive variables are taken into account. For 
the entropy O we have 



O = -X-yj^obsX 



= -^ + ^A^.] (40) 



dX, 



0{y,} = -X{y,)~y,x, (41) 

For the case of a completely open system one has X = — O and, hence, from Ea. (|36|l 

O = In ^ g{x.}e"«'^' {2/0 = {yi} (42) 



and so, in general, for any system we have 



O = In ^ 5{x,}e-^'^' (43) 

{X.} 



In the case of macroscopic systems, O = O^y.y = 0. All above equations can be supplemented to the following, relating 
the number of characteristic elements of a system to the countable microstates resulting from specifying extensive 
variables. 

9{x,}u{y.} = J2 5{x,}u{x.}e-^'^^ (44) 

{X,} 

We have deal with discrete specifications of the elements within a class. If a continuous spectrum is considered, the 
last equation must be written as 

r -r-r dX 

9{x,}u{y.}= ff{x.}u{x.}e-«-^' J]— ^ (45) 

JiXi} j '-^Xz 

We see that in the case of a continuous spectrum the Laplace transform is the appropriate mathematical tool for 
passing from one ensemble to another with less extensive constraints. The reverse process can then be mathematically 
described in terms of the inverse Laplace transform. By using the Bromwich inversion formula [l9l | , we can write 

9{x,}u{x,} ^ IT- 5{x,}u{i;.}e^'^'n7r^ (46) 

where 7{j,i} is chosen so that the poles and singularities of the integrand in the complex plane are at the left of the 
integration path for each of the variables {yi}. For small systems, however, these transforms become non-local leading 
to ensemble non-equivalence. Eas. (|36|) . H37|) . H38|l . (|39|l and H43|) are the most important expressions of Boltzmann- 
Gibbs equilibrium thermostatistics. 
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D. Probability 



From the previous treatment we can define now the probability V of measuring values for extensive non-environment 
variables {Xg} {{Xg} C {^i}) as the chance of the system for being in a given configuration compatible with the 
specification of these variables. This is simply the quotient of the numerical value of configurations at a specified 
value of variables {^s}, g{Xj}\j{yi]'j{x^} to the value of the whole collection 



V. 



{Xs} 



9{X,}u{y,]U{X,} _ 9{X,}u{y,}U{X,} 



9{X,}V{y,} E{X,}ff{X,}U{y,}U{XJ 

,9{x,}u{x,}e-f=^= ^^^^ 



T.{Xs}9{x,}u{x.]e y^^^ 



From Eas. l|^ and we obtain 



V{x,} - e^-^t-=> (48) 
The probability distribution is normalized, as is evident from Ea. l47|) . 

^ e^-X{x,) = I (49) 

{Xs} 

The values of non-environment variables fluctuate as these are not specified a priori. The observed macroscopic 
values for non-environment variables are given by the differential equation H37|l obtained from thermodynamical 
considerations. From that equation, we have, for a variable belonging to the set {^i}, 



Xs,obs = 7^ (50) 



{X,},{y,}i^s 



From Ea. (|5n|) and using Eas. (|^ and and the definition of probability given above, Ea. l|T7|l . a straightforward 
calculus leads to 

Xs^obs — — — < Xs > (oij 

i.e. the observed value of the variable Xg is the ensemble average of this quantity < Xg > over all attainable 
configurations. If the observed values are related to time averages, Eq. H51|l constitutes the so-called ergodic hypothesis 
which is of fundamental importance for the probabilistic approach in Boltzmann equilibrium statistical mechanics. 

The probability Pk.{Xs} for the system to be in a concrete configuration k compatible with the specification of the 
constraints and variables {Xg} {{Xg} C {Xi}) can now be defined as 

T^{Xs} 9{X,}u{y,}u{Xs} 
Pk.tXA = = 

9{x,}u{Xs} 9{x,}u{Xs}9{xi}u{y-} 
e-VsXs 



E{Xs}9{x,}u{Xs}e y^^' 

p-VsXs 

(52) 



'VsXs 



In terms of this probability distribution, the dimensionless characteristic thermodynamic function, Ea. (|36|) . can be 
rewritten as 

A- = -ln^e-^'-^'l'= (53) 

k 

From Eqs.ljSOI) and we obtain 

Pk.lXs} = e^-^t-=>--^t-=> (54) 



11 



which is also normaHzed leading trivially to the same result concerning observable values, Ea. H51() . For the entropy, 
we have, by using Eas. (|^ and 

J = ViXi^obs - X = '^pk iviXi - X)i^ = 

k 

= T.p'^^^^S^ = -T.p'^^^p^ (55) 

k ' fe 

which is the Gibbs-Shannon entropy. It has been used by many authors, since the pioneering work by Jaynes (3], to 
develop ensemble theory. 



E. Theory of Fluctuations 

In this section we follow an analogous treatment to the proposed by Miinster to develop fluctuation theory [20l| . 
We consider macroscopic systems in this section (and, hence, = 0). Ea. (|^ allows to calculate the probability of 
having a concrete value for a fluctuating set of macroscopic parameters {a}. To simplify notation we define the vector 
a = (ai, q;2, ...) containing the fluctuating variables. We have, from Eq.|0H)l 

■Ps - Vsoe'^^o-^^ (56) 

Where do is the most probable value of the variable a at an equilibrium state. We take q?o = for simplicity (but 
with no loss of generality). When fluctuations are small, a Taylor expansion of A3 in Ea. (|55|l to second order yields 



Vs = CPoe ^y^-^-'Jo (57) 

where C is a normalization constant and subindex denotes now the final equilibrium state. Ea. (|57|l is the so-called 
Einstein fluctuation formula. It can be applied to calculate the variances of fluctuating macroscopic parameters. In 
Ea. (|56|l AX = — [Xs — X{0)] is to be considered as the (dimensionless) work W developed by a subsystem on its 
surroundings in evolving to an equilibrium state 0, IT^ . This work is developed by keeping intensive environment 
variables constant, and, hence, from Ea. H37|l 

AX ^SW ^ -XAd (58) 

where A must be considered as the vector of parameters conjugate to the vector a. Note that Ea. l|58() makes reference 
to a reversible process, in fact W can also be pictured as the minimum work Wmin which is necessary to do on the 
subsystem for driving it from the equilibrium state to the nonequilibrium one. As this work is the minimum one it is 
developed by following a reversible path It is useful to define the matrix G as 

1 

Gki^(^] =-(,^1 (59) 



\dakdaijr. \ dXkdX. 




where in the last equality X^ is defined through a Legendre transformation of X by changing variables ak (or ai) to 
Afe (or Xi) in Ea.(|S5|l 

AX+ = -X*dd* + audXk (60) 

The asterisk in Ea. (|60|l denotes vectors with vanishing k position. From the above considerations, Ea. (|57|l can be 
rewritten in the following way 

and, therefore, from Ea. fK^ and Eq.lj^ 

A = Ga = -^ (62) 
oa 
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For the average of the dyadic product aX (in components {dX)ki = OffcA/) 

< dX >= / dXVadct = d dd (63) 

and by means of integration by parts and and considering the normahzation condition for the probabihty distribution, 
we obtain 

< aX >= U (64) 

whence 

< AA > = G (65) 
<dd> = (66) 

The following important equality holds 

< AA >< dd >= 1 (67) 
If a is chosen to represent the extensive variable Xs , we have, from Eq. (|66(l and by using Eq. (|37|l 



dX^ Jo V% 

for the mean square fluctuation. If the conjugate intensive variable is now considered we obtain, from Eqs. H65|) and 
(|57|l the following result 

<-">-(t)>-ffl: 

where superindex + has been dropped in the last equality for X which must be understood as the concrete specification 
of this function for the set of variables chosen (those involved in the derivation plus those held constant). 
The Ea. (|67(l for the extensive- intensive variables pair reads 

< AX^ >< Ay2 >= 1 (70) 

when the variances are evaluated by keeping constant the same environment variables. If there are several extensive 
variables that can fluctuate the covariances among them can be easily calculated from Ea. (|57|l . The result for the 
mean square covariances of two variables ^'^si, is 

<AX.,AX,,>= U—^i-^ (71) 



dXs.dX, 



S2 







F. Expressions for the customary statistical ensembles. 

In this section we apply the main expressions obtained previously to the customary statistical ensembles. It is useful 
to abbreviate notation to define the following quantities 13]: (3 = tt = v = The following expressions 

were also obtained in the previous formulation by Vives and Planes ]13| . 

1. Microcanonical Ensemble 

For this ensemble, {X^} = {E,V,N}, {y^} = {(3,tt,v}, {Xi] = {0} and {y,} = {0}. For this system Eqs.ljSll) and 
()38|l take the form 

X = ~^ = -\nn (72) 
k 



13 



The values for non-environment variables can be calculated from Ea. (|37|l 



(3 



dX 
'dE 

dX 
dV 

dX 

9iV 



V,N 



E,N 



E.V 



The probability distribution for a configuration pk is given by Eq. H52|) . 



The variances of non-environment variables are given by Eas. (|68|l and (|69|l 



EM 



df3_ 

2 v\ -1 



d^X 



E.V 



d^xy' 



(73) 
(74) 
(75) 

(76) 

(77) 
(78) 
(79) 



Note, that on evaluating the fluctuations, the independent variables considered are those held constant and the one 
involved in the derivation. X is considered as a function of these independent variables and hence, Ea. (|37|l is used 
each time with different sets of independent variables. For example, on evaluating the variance of variables (3, V, N 
are taken as independent in writing Eq.(j2Zl, but, on evaluating the variance of tt, the independent variables are tt, 
N. This means that X is different in each case, but if one uses the general definition Ea. H37(l for X, these differences 
are of no formal interest as a complete set of environment variables specify those which are not natural in the system, 
through Ea. (P7jl . in general. 



2. Canonical Ensemble 



In this case {Xj} = {V,N}, {y.j} = {tt, z^}, {X,} = {E} and {yj = {/?}. For this system Eqs.(|21l and (EHl take 
the form 

X ^ PF =^ (3{< E > -TS) = -InZ (80) 

where F is the Helmholtz free energy and Z = g{f], N, V) is the partition function of the system. 
The values for non-environment variables can be calculated again from Ea. (|37|l 

'dX 
dV 



<E> = { — ] (81) 

V,N 

(82) 
(83) 



dX \ 

m) 



The probability distribution for a configuration Pk,E in which the system has energy E, is given by Ea. H52|l which in 
this case takes the form 



Pk,E = (84) 
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The variances of non-environment variables are given by Eas. (|68|l and (|69|l 



dE 



,2 --^^ -1 



(85) 



(87) 



3. Grand Canonical Ensemble 

For a system that can exchange energy and particles with a heat bath at T and fi, we have {Xj} ~ {V^}, {%} = {^}: 
{Xi} = {i?, A^} and {y^} = {/3, i^}. For this system Eqs.ljSEI) and (j2Hl take the form 

X ^li{<E> -TS) + v <N >^-\i^-~ (88) 

where S = (7(/3, V) is the grand partition function of the system. 

The values for non-environment variables can be calculated again from Ea. (|37|l 



<E> = \ ^-^] (89) 

V.y 



dp 



(90) 



<^>-(f)^ (91) 



0-y 



The probability distribution for a configuration Pk,E,N in which the system has energy E and number of particles N 
is, from Ea. l(5^ 



^-pE-vN 



Pk,E,N = = (92) 

The variances of non-environment variables are given by Eas. (|68|l and (|69(l 



(95) 



Since for this ensemble there are two extensive fluctuating variables, we can calculate also the correlation between 
them, from Ea. l(7T|) 



< Aii;AiV> ^ f^^V' 
\dEdN J 



V 

dN \ _^_fdE\ 



(96) 



III. SQUEEZED STATISTICAL MECHANICS: A GENERALIZATION OF THERMOSTATISTICS TO 

SYSTEMS OUT OF BG EQUILIBRIUM. 



Hitherto we have considered systems with no capability to choose between the attainable configurations compatible 
with the constraints. The systems treated here fill their attainable configurations depending on their degree of 
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attainability, as specified by the physical constraints imposed to the system. The features of the constraints have 
been the only criterion for the system to decide what configurations adopt. As we have seen these constraints can be 
extensive or intensive, and are imposed to the system thus restricting its behavior. 

There are complex systems in nature that exhibits self-organization and have the ability of taking decisions. These 
systems can restrict (or enlarge) the total number of configurations at their disposal by selecting only a group of 
them or expanding their spectrum against constraints. An analogous case can be that of systems in nonequilibrium 
situations for which the total number of configurations of the phase space are not necessarily available. 

We thus see that new information is required in order to account for these systems. A possible way to introduce 
this information is to modify the definition of X by introducing not the complete set of configurations but the actual 
one g = hg(g). The function h specifies the information of the actual set of attainable configurations for the system 
in reference to those attainable from the only specification of the constraints (g) and then act as an interface between 
equilibrium and complex systems or nonequilibrium situations. The parameter q is introduced for the sake of generality 
to classify different statistics under the common denominator of a superstatistics, specified by the squeezing function h. 
We remark that the function h can have a more complex dependence on system variables and situations. We consider 
here, however, the case of one-parameter dependent squeezing functions, in which the complexity is accounted for 
through an entropic parameter q and the specification of the function h. Obvious extensions are considering n- 
parameter dependent squeezing functions hq-^,,,,q^. A theory for these can be developed following analogous lines 
to those we discuss below. The convenience of introducing at least one additional entropic parameter stems in the 
diversity of the dynamics of complex systems. We introduce the function kg as an additional postulate 

Postulate 3 The actual classes of configurations g of a complex system can be determined by a convenient positive 
squeezing function hg{g) of the BG classes g by means of the relation g = hq{g). This function has a well defined 
inverse function Hg, being hq{Hq{x)) — Hg(hq(x)) = x. 

The set of fundamental equations H3()|l . (|37|l and (|38l) can be rewritten then as 

X ^ -lng{x,}u{y.} = -ln/ig(g{x,}u{y,}) (97) 
dX = -Vj.obsdXj + Xi^obsdui (98) 

J = -X + UiXi^obs = -X + yi [ —— 



ln5{Jf,} = ln/ig(g{x,}) (99) 

dX 



O ^ -X ~y,^obsXj = -X + X 



9^3/ {x,,},{y^} 



= ^'^9{yi] =^'^hg{9{vi]) (100) 
9 = hq{g) (101) 

If we compare these expressions with Eas. (l36|l . (|37ll . (|38|l and H4U|I . we see that the modification introduced concerns 
only the attainable configurations by the system through Eq. Hl()l|l and anything else. Additive decomposition of 
the actual (squeezed) classes (as the fulfilled by the Boltzmann-Gibbs equilibrium classes) does no longer hold in 
principle, however. Expressions like Eas. H19|l and H20|l cannot be written for the "squeezed" classes 'g. This fact has 
its deep grounds in the underlying nonlinear kinetics exhibited by these systems (see below). 

We call the framework described by Eqs. (|97|l to (|101|l Squeezed Statistical Mechanics. The term "squeezed" 
makes reference to the function hg that deforms the Boltzmann-Gibbs equilibrium phase space. The expressions 
above are valid for any statistical system and can be considered as the skeleton of the thermostatistics presented 
here. Nonextensive thermodynamics demands O 7^ but it does not necessarily imply a nonextensive statistics. 
Nonextensive statistics is obtained when the function hq is different to the identity. In the case of extensive statistics 
Eqs. (jnZJ) to H1U1|I are equal to Eas. (|S^ . (|S7|l . and (gUl with hq{x) = x. The main motivation to introduce 

hq{x) ^ a; is a practical one: it is usually more feasible counting over the complete distribution of configurations 
(i.e. the equilibrium Boltzmann closure) than over the incomplete distribution which constitutes the nonequilibrium 
closure. When this is done, Boltzmann's principle cannot be applied with the equilibrium closure alone and an 
additional information is required to relate the actual total number of configurations to that of the equilibrium closure. 
The different statistics arise when one uses the equilibrium (Boltzmann) closure to establish probability distributions for 
complex systems or nonequilibrium situations. As a result of that, the equilibrium closure and classes of configurations 
are squeezed and rearranged, and no longer resemble those of BG equilibrium. We show this point below. 

Gibbs fundamental equation written in the generalized form, Eq. H98|l , is formally identical to the obtained above for 
the ordinary thermostatistics, Ea. (|37|) . The robustness of the Legendre transform structure and Gibbs fundamental 
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equation for arbitrary entropic forms and mean values has been recently pointed out in the references |22| in the 
framework of the Jaynes approach. 

The actual total number of configurations Ij is then to be understood as dependent on the total number of configura- 
tions g compatible with the constraints and, maybe, on new parameters. We emphasize the fact that configurations Ij 
represent, in any case, the actual configurations of the system and not those resulting from counting over the complete 
distribution g which are related to a system under the hypothesis of the molecular chaos (i.e. with no organization 
and no ability of taking decisions) at equilibrium. This case has been considered in previous sections. Again, the 
methods used in previous sections to formulate equilibrium statistical mechanics can be identically translated here to 
the extension to nonequilibrium situations and complex systems. 



A. Generalized ensemble theory. 

Let us write the squeezing function in terms of the actual and the BG characteristic classes 

9{X,}vj{y,} = hq (g{x,}u{y.}) 
The function hq is equal to the identity for BG statistics. We have, from Ea. (|97(l 



X - -\nhq {9{x,}vj{y.}) 



and thus, if Ea. l99|l is used 

^{x,} = ViXi - J{x,} = ViXi - \nhq {g{x,}ij{Xi}) 

hence, from Eas. HlU4|l and (|105|l 

hq (5{x,}u{y.}u{x,}) = hq {g{x,}^J{x,}) eT^'^' 
Let us consider the inverse function of hq^ Hq, so that Hq{hq{x)) — x. We have 

5{x,}u{y.}u{x.} = Hq {hq {g{x,}u{x,}) e~y^^^) 

and, therefore. 



X 



In h„ 



■ViXi 



(102) 

(103) 

(104) 

(105) 
(106) 
(107) 

(108) 



This is the central equation of generalized ensemble theory. It has to be supplemented to the following for the two 
entropic forms that we have been considering J and O, the last being of relevance for small systems. 



J = ln^9(g{jf^}) = \nhq 



Yl 9{x,}yj{x, 



O = ln/i,(g{,y.}) = \nhq 

All these expressions equal to those of equilibrium when the squeezing function is equal to the identity. 



-v.x. 



(109) 
(110) 



B. Extensive and nonextensive representations. Subdivision property 



The above results suggest that there exist two convenient representations when nonequilibrium or complex systems 
are considered. These two representations are equal only for systems following BG statistics. We have, on one hand. 
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the extensive representation specified by the actual classes g. In terms of these the entropy, specified by Eq. H99|l . is 
extensive 

Ja+b hi5A+B = IngAffB = Ing^ + Ings ^ Ja + Jb (111) 

From this extensivity property we have another very interesting property. The subdivision of the actual classes of a 
composite system in BG classes of independent subsystems has the following form 

9A+B = hq{gA)hq{gB) (112) 

We call this property of factorization an actual class in BG subclasses subdivision property of the actual class. 

We have also a representation specified by the BG classes which leads to nonextensivity in the entropic form (as 
BG classes are not the proper representation of the system). This can be clearly seen from Eq. H99|) . If the inverse 
function Hq is used one straightforwardly obtains 

Ja+b = InhqigAgB) = Inhq [Hq{e^^)Hq{e^^)] ^Ja + Jb (113) 

Although the BG classes are not the proper classes of the system, it is very convenient working with them due to the 
very useful properties these have and which are lacking in general in the proper actual classes, the property of linear 
decomposability Ea. (|20l) being an example. Furthermore, the equilibrium classes can be known better as they belong 
to the complete description of the BG closure. 



C. Probability 

We maintain unchanged the definition of probability we given at equilibrium (see section III Dl) . The reason is 
that we usually use the BG classes to evaluate probabilities since these constitute the complete closure of the phase 
space of the system. The form of the probability distributions is no longer that of BG thermostastics, however: the 
classes are rearranged as we show below giving rise to bunches of distributions that can be classified within a concrete 
superstatistics. Each of the distributions has to be understood as a different way of partitioning the same equilibrium 
closure. For the probability we have, then 

„ 5{Xj}u{y,}U{X4 5{X,}u{y.}U{JfJ 



5{Xj}U{y,} '^{X,} 5{Xj}u{y.}U{X4 

Hq{hq (.9{x,}u{x.})e-^=^0 



T.{x^} Hq {hq {g{x,}ij{x^}) e-«-^=) 



(114) 



H,('i,(s{x,}u{x,|.})e-''=^=l'=) 



Pk.{x,} — - — = T-f — N _„ X (115) 

9{Xj-}u{Xs|fc} 

where Xs\k denote the value of variable Xs in the configuration k and g{Xj}u{Xs\k} the total number of microcanon- 
ical configurations within the class to which the configuration k belongs. The generalized Boltzmann factor B^Xi} 
can then be defined as 

_ 9{X,}U{y,}U{X,} Hq (hq }U{JC, }) g-^-^') 

B{Xi} = = (lib) 

9{Xj}u{x,} 9{Xj}u{Xi} 



From Eas. HlU5|l . (|108|l and (|109|l we obtain 



i/,(e-^f-=}) 



Hq{e--^i^=\t'}) 

^'^■^^^^ = Hq{e^i^s^^y)Hqie-^) ^''^^ 



These probability distributions are normalized. 
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For the observed values of extensive non-environment variables, we have from Eas. H98|l . (|108|l and H115|) and after 
a straightforward calculation that 

V%/{x,},{a,.} /(^{x>}5{x,}u{ya) 
/(g{x,}u{jy.}) Y 

^ g{x,}vj{x,^^}}{Pkg{x,]u{x,\^}g{x,}vj{y,]) 

where /(x) = dhq{x)/dx is the first derivative of the squeezing function. We note that the definition of the mean field 
value is very different from that of BG equilibrium although the same definition for probability has been used. This 
is due to the fact that the classes of configurations of the equilibrium closure have been rearranged and the partition 
is different now. 



D. Connection with Beck-Cohen superstatistics 

A very recent approach to generalize statistical mechanics to nonequilibrium from the bottom up is that of Beck 
and Cohen |25||. The following discussion will be referred to the canonical ensemble only, as is done in Beck-Cohen's 
paper. The approach followed by the authors is based on the definition of the generalized Boltzmann factor B{E) 

poo 

Be = / d/3'f{P')e-^''' (120) 

where /(/3') is a function of the thermal Lagrange parameter /?' which is a normalized probability density that allows 
to build a normalizable statistics which must be reduced to BG-statistics when there are no fiuctuations in /3'. In 
fact the choice /(/?') — S{(3' — (3) allows to recover the ordinary Boltzmann factor e~^^ . We can relate easily the 
formalism of Beck and Cohen to the presented here. From Eas. (|116f) and H12UI) we obtain 

/(/?') = L'^' 

where L^^ denotes inverse Laplace transform with respect to variable /3. We see that an interesting property (but 
not, in fact, a requirement) for the function H is to have a well defined inverse Laplace transform. 



Hq {hq {g{N,v}u{E}) e 

g{N,V}u{E} 



(121) 



E. Zeroth law of thermodynamics: a pathway to the relevant statistics. 

Within the extensive framework, if two systems I and II are put in contact attaining thcrmodynamical equilibrium 
their Lagrange parameters become equal 

- yf (122) 

From Eas. (|97|l and (|98|) . this equality can be rewritten within the nonextensive representation (microcanonical en- 
semble) as 

rfing^ ^^l.BG = y^BG figx, )yiBG = figx, )yj!BG (123) 

where yj,BG = ~dhig^Xj}/dXj is a BG Lagrange parameter which, if the thermal one is considered, can be measured 
by an ordinary thermometer following BG thermostatistics. We arrive to the following conclusion: the temperature 
measured by means of a BG thermometer for systems following the generalized zeroth law principle is, in general, 
different for each of the systems in contact if those are of different nature (i.e. being described by different statistics). 
This is not surprising since temperature measured by a thermometer depends on calibration, having a statistical origin, 
and we are describing non conventional situations out of BG equilibrium hut following a more general definition of 
thermal equilibrium in which only temperatures measured by proper thermometers are equal. If a system does not follow 
BG statistics, a proper thermometer should be developed to measure the actual temperature at equilibrium. Stated 
in other way, for two systems in (general) thermal equilibrium, the proper thermometer of each system measures the 
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same temperature but temperatures measured by BG thermometers are different (in general) if those are of different 
statistical nature. 

Let us consider now that the system denoted by II is described by BG statistics {f{gx. ) = 1) and it is in equihbrium 
with system I which is a complex one. Eq. (|123|l can be integrated to give 

hg{9'x,)=K{g'x])e""''^ "^'^^ ' (124) 

where hq{gx) is a reference state. If we take hq{l) — 1 when there is only one attainable microstate, having no 
squeezing in this case (which seems reasonable as all superstatistics should then collapse) and we take series of 
measurements of yj,BG for both systems in (generalized) thermal equilibrium calculating their ratio as a function of 
the number of BG attainable microstates of the complex system, Eq. (| 124(1 provides a mechanism to determine a priori 
the statistics that the latter follows. 



F. Irreversibility and generalized Boltzmann equation 

The actual attainable classes of microstates g introduced here can be directly related to the underlying nonlinear 
kinetics in the squeezed statistics under consideration . The kinetic interaction principle first stated by Kaniadakis 
p?! leads to a generalization of the Boltzmann's kinetic equation allowing to consider more complex systems. Our 
formalism is completely consistent with the previous Kaniadakis theory, which in fact, sheds some light on the meaning 
of the squeezing function hq. 

The number of particles in a volume drdv (where r and v are the position and the velocity of a given position inside 
a fluid) is given, at time t in terms of the one-body density function F{f,v,t), by the product F(f,v,t)drdv. If we 
use the above established subdivision property for the density F we can assume that the number of pairs of particles 
in the configuration space element dr and velocities dvi and dv2 is given by the following relation 

hq[F{r, Vut)]hq[F{r, V2,t)]dvidv2drdf (125) 

This assumption is reasonable if we realize that v specifies a general actual microstate of one system and we can thought 
this system to be subdivided in parts at different microstates counted in the BG closure with velocities vi and 1)2 ■ 
In the case of BG systems (where hq is the identity) the assumption coincides with the celebrated Stosszahlansatz 
(molecular chaos hypothesis) due to Boltzmann. It cannot be justified from first principles and it is certainly nor 
necessary nor true for all time |28j| . We are considering here a more general form for the number of pairs of particles 
arising from squeezing of the different elements of phase space and leading to a general factorization of the one particle 
distribution function. It also leads to a very general nonlinear interaction which has been discussed previously I23 
and from which the following generalized Boltzmann equation is obtained (if we replace k(/) in the discussion of |27| 
by hq{F)) 

d 1 dV{x) d 



dt 



^F', F)^{F[,Fi)T{t, X, V, V', iTi, «])} (126) 

where V{x) is an external potential, x is the position vector of a site in the fluid and, v,v\vi,v'i are the initial 
and final velocities of particles on this site and on another one being both in interaction. F, F', Fi and F[ are the 

initial and final particle densities of each site. T{t, x, v, v' , v[) is a factor which depends only on the nature of the 
two body particle interaction (proportional to the cross section) and ^{F,F') = £^{F' , F) is an arbitrary symmetric 
function. This equation reduces to the traditional Boltzmann equation if hq{F) = F and ^(F', F) — ^{F[, Fi) = 1. 



d 1 dV{x) d 



V 



dt dx m dx dv 



F = / d''v'd''vid''v[[F'F[~ FFi]T{t,x,v,v',v{,v'^) (127) 



Importantly, if hq{g) is the identity, the variational principle leads to the Maxwell-Boltzmann distribution regardless 
of the value of the function ^. This result indicates that there is an infinity of ways (one for any choice of ^) to obtain 
the Maxwell-Boltzmann distribution. 

If ^ is the chemical potential of the particles and U{v) their kinetic energy, the generalized Boltzmann equation 
implies that the Lyapunov functional Ti given by [23 

H = - / d"a;d"?; / dF\nhq{F) - /3 / d"a;d"w [Vix) + U{v) - fi]F (128) 
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is non decreasing with time if f{F) = dhq{F) / dF > and hence as the second integral in the r.h.s of this equation 
is a conserved quantity, the entropy S = — J^d"xd"v J dF\iihq{F) satisfies also that dS/dt > 0. The condition 
f{F) was advanced implicitly (for positive BG temperatures) when studying the generalized zeroth law (see above) 
provided that hq is a positive function in Eq. (|123|) . We remark that being hq an arbitrary function, H-thcorem is thus 
verified for a very large class of nonlinear systems. This kinetical picture completes the above formulation of Squeezed 
Statistical Mechanics, in which a general thermodynamic framework for extensive/nonextensive macroscopic/small 
systems has been provided. 

G. "Squeezing and rearranging" 

In the above sections we have formulated BG statistical mechanics through a new method which can be called 
"rearranging". Related to it are the concepts of BG classes of configurations, linear decomposability and of the 
existence of a function for any system which is a monotonic decreasing function when increasing elements in the 
characteristic class. The key step in the method is that which goes from Eq. to Ea. H33|l . 

In extending BG thermostatistics to other many relevant situations we have shown that a previous step is needed in 
order to define the statistics. This previous step is called "squeezing" and consists on the specification of the function 
hq and calculating (if possible) its inverse Hq and its first derivative /. Then a general version of the rearranging 
step, Eqs. H104|l to H106|l leads to the decomposition of the equilibrium classes and, hence, to a workable version of 
the non conventional statistics at hand. 

IV. APPLICATION: TSALLIS THERMOSTATISTICS 

A recent important effort to extend statistical mechanics to complex systems has been due to Tsallis . From ideas 
taken, in part, from multifractal analysis and information theory, Tsallis proposed a definition for the entropy based 
on the parameter q and after applying the condition of extremum for the entropy, he provided an extended formulation 
of statistical mechanics. In the following paragraphs we develop Tsallis statistics on the context of our framework 
by introducing a particular choice on the q-squeczing function for the total attainable configurations compatible with 
the constraints. 

Let us consider a complex system (I) in (generalized) thermal equilibrium with a BG thermal bath (II) at Pbg.ii- 
Measurements with a BG thermometer of Pbgj are taken over the complex system by varying (for example) its 
volume. Let be the BG microstates of the complex system known each time. These are all those allowed by the 
constraints, that can be controlled externally. Let us further consider that a power law (we abbreviate by gj the 
number of attainable microstates of the complex system) controls the quotient of Lagrange parameters (which is not 
unreasonable due to the ubiquity of power laws in nature) 



Pbgj 
Pbgji 



the exponent 1 — g is introduced here for convenience. If q = 1 we then would have BG equilibrium between two 
BG systems. If this power law is introduced in Eq. 1124() and we integrate from hq{g]^,) = hq{l) = 1 we obtain the 
following form for the squeezing function hq{x) 

/ig(x)=e''T^ (130) 

The inverse function Hq{x) and the first derivative fq{x) can be computed straightforwardly 

Hq{x) ^ [(l-q)lna; + l]i/(i-9) (131) 
fq{x) = x-^ (132) 

This completes the "squeezing" step. The knowledge of these three functions is enough to develop the thermostatistical 
framework. We use them in all what follows. The following relations hold 

ln{hq{x)) = lliqX (133) 

Hqien = e^q (134) 
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In these expressions In^ x and are, respectively the q-logarithm and the q-exponential functions first defined by 
Tsallis 

a;^~« - 1 

ln„ X = (135) 

1 - q 

e^„ = [(1 -<7)x + l]i/(i-9) (136) 



A. Entropy, characteristic function and ensemble theory. Nonextensivity 
From Eqs. ()108|l and (|1(J9|I we immediately obtain 

X = JjSmvAZI (137) 

J = (138) 

1-9 

for the characteristic function and the entropy, respectively. The form of the characteristic class g{Xj}u{yi} can be 
obtained easily from Eq. H1U7|I 

E'^-^ / yiXi{q — 1) \ ' 
.9{x,}u{j/.}u{x.} = 2^ 5{Jf,}u{x.} 1 + (139) 

{Xi} {X.} \ 5{x^.}u{x,} / 

We observe that, for the microcanonical ensemble, Tsallis microcanonical partition function coincides with the 
Boltzmann one (fi), as it is easily seen from Ea. H139| by considering that for a system in the microcanonical ensemble 
{X,} = {y.} = W: 

The nonextensivity property for the entropy in terms of the BG classes can be also straightforwardly obtained from 

Eq. (ins 

Ja+b = \nhq{gA9B) = In/ig [Hq{e^^)Hq{e^^)\ ^Ja + Jb + [1- q)JAJB (140) 



B. Probability 

The probability distributions can be obtained from Eqs. 1114(1 and ((115(1 



l-<! 



V{x.} = ^ lif,}™/ ^=e^'1 J (141) 

' ^ 1-9 



E{x,} .9{x,}u{x,} ( 1 + 



1 , y.x,U(g-i) 



1 + ^{X3}(9-1) 



Efei 



ysXs\k(q-l) 
T^^i 



where Eqs. I(117() and 1(1 181 have also been used. These probability distributions are normalized. 
The observed values of extensive non-environment variables can be obtained from Eq. ((119(1 



Pi 



X, k [9{Xj}\j{x,\k}\ 



JTr^Xi]]^ — < Xi >q 



where < Xi >q denotes the q-mean value. In general, for a mechanical quantity A, we have 



(142) 



(143) 



<A> 



1 \ li-g 

fc [9{Xj}yj{Xi\k} 



{X,}yj{yi}lS{X,\k} 



(144) 
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This q-mean value is different to those suggested by Tsalhs et al. |2J] . Let us analyze briefly some of the implications 
of this result. In particular, we see that in the microcanonical ensemble the norm is conserved, i.e. 



< 1 >a 



E 



9U} 



9{Xi}9ix,} 



(145) 



The norm is not conserved, in the rest of statistical ensembles. Only integrals of motion are properly defined in the 
microcanonical ensemble, however, and the lack of normalization in other ensembles is not a matter for worrying 
about here. The average J2k Pk^\k / J2k Pk recently considered could be obtained with other definition of probability 
different to the given above. We choose however to maintain the same definition of probability as we given for 
equilibrium. Operationally, the mean value defined by us has the same advantages that the so-called second choice 
for the mean value proposed by Tsallis: < A >q= X^fc-Pfc^U- Furthermore, the conservation of the norm is at least 
preserved in the microcanonical ensemble. 

For the entropy, we have also have, by using Eas. (|^ . H137|l and 1] 143(1 



J7 — yiXi,obs ~ ?i 



1 I yiXi\k(\-q) 

^^^ \ k [5{x,}u{x,U}] 



-1 + 5{4}u{a.} ^ 



Pk ?'fc-9{x-}u{y.}J +5'{x^}u{y.} 1 



1-9 



(146) 



which is the Tsallis generalized entropy. This was the starting point for Tsallis to develop nonextensive statistics. 

Einstein's fluctuation formula can be easily deduced for Tsallis statistics by considering a general vector a in 
Ea. (|141f) instead of variables {Xs\ and expanding to second order in a 



In 



We obtain 



(147) 



(148) 



If the matrix G defined above, Ea. ((59|l . is used, the variances can now be calculated from Eg. 1(148(1 for a. and A = Ga 
following an analogous procedure to that of Section III El These are given by 

(149) 
(150) 



<AA> = [\ + {q-l)Xo\G 
<aa> = [l + iq- l)Xo] 



And thus 



V< AA >< da > 



(151) 



It is noteworthy that in the limit g ^ 1 all expressions obtained in this Section reduce to those obtained in Sections 
III CI III Dl and lTlEl It must also be noted that the thermodynamic formalism exposed in Section III Al does not change. 
A clear connection between Tsallis statistics and thermodynamics is thus established and Legendre transform is 
warranted in the thermodynamic limit. 

The treatment followed embodies Hill thermodynamics and, hence, is stable, allowing for large fluctuations in 
physical quantities. The formulation has been based in the concepts of classes of configurations and in the principle 
of minimum for the dimensionless (generalized) characteristic thermodynamic potential. Both concepts have a clear 
meaning in thermostatistics as we have shown. 

We conclude remarking that we have formulated a possible generalization of BG thermostatistics in Section ITTll 
We call this formalism Squeezed Statistical Mechanics due to the fact that it makes reference to a squeezed phase 
space because of the introduction of the squeezing function hq. In the framework of this alternative formulation we 
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have developed classical equilibrium, BG statistical mechanics (supposing no squeezing) and a general framework 
valid for nonequilibrium situations and complex systems. Among the advantages of the formulation presented here 
are the compactness and generality of the expressions obtained and the rationalization of recent statistical mechanics 
developments in an unified approach. New insights in these developments have been provided. Furthermore, the 
general formulation presented is fully consistent with a generalized kinetic theory previously proposed. Applications 
to concrete physical systems of the formalism presented will be given in forthcoming papers. We indicate however that 
previous applications of Tsallis and Beck-Cohen statistics and also Hill's nanothermodynamics can be considered also 
as applications of the extended theory presented here. Moreover, the treatment proposed can be extended to other 
more complex choices for the squeezing function different to the one-parameter dependent family we have considered 
in this paper. Linkings to the very interesting and important formalisms including SRB distributions and the chaotic 
hypothesis are also to be explored. 
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for several comments on the manuscript. One of us (V. Garcia-Morales) wishes also to thank support from the Spanish 
M. E. C. D (grant No. AP2001-3329) and conversations with J. M. Garcfa-Sanchis. 
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